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ABSTRACT 

The  mode  conversion  of  an  electrostatic  whistler  wave  into  a 
Bohm-Gross  mode  at  plasma  resonance  is  analyzed  fo^  a  magnetized  plasma 
with  a  longitudinal  density  gradient  (i.e.,  Vn^  x  B  -  0) .  It  is  found 
that  a  whistler  incident  upon  plasma  resonance  from  inside  the  plasma 
converts,  without  producing  a  reflected  wave,  into  a  short  wavelength 
Bohm-Gross  mode  that  carries  energy  down  the  density  gradient  away 
from  resonance.  The  detailed  structure  of  the  electric  field  near  the 
resonance  is  found  analytically.  It  is  shown  that  the  production  of 
the  Bohm-Gross  wave  by  mode  conversion  can  be  described  by  a  model  of 
plasma  resonance  driven  by  a  k-0  electric  field  (i.e.,  the  capacitor 
plate  model).  The  relation  between  the  driver  amplitude  and  the 
amplitude  of  the  incident  whistler  is  derived. 
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I .  Introduction: 

This  paper  analyzes  the  process  of  mode  conversion  of  a  long 
wavelength  electrostatic  whistler  wave  into  a  short  wavelength  Bohm-Gross 
mode.  The  whistler  with  frequency  w,  less  than  the  electron 
gyrofrequency ,  Q  (i.e.,  w  <  n) ,  propagates  from  inside  the  plasma  in  the 
direction  of  decreasing  density  to  the  plasma  resonance  point,  Wpg  -  w, 
(wpg  is  the  electron  plasma  frequency)  where  mode  conversion  occurs. 

The  electrostatic  whistler  waves  considered  here  satisfy  the 
dispersion  relation  _  q,  where  in  the  cold  plasma 

limit  £||  =  1  -  (jpg^/w^  and  -  1  -  .  In  a  magnetized  plasma 

in  which  the  density  varies  only  along  the  magnetic  field  direction,  z,  as 
illustrated  in  Fig.  1,  the  WKB  wave  number  component  along  the  magnetic 
field,  k||=  k^(-£^/£:||)  ,  is  a  function  of  position.  As  the  whistler  wave 

propagates  towards  lower  densities  so  that  Wpg  approaches  u,  £||  approaches 
zero  and  the  parallel  wave  number  becomes  large.  The  inclusion  of  finite 
temperature  effects  in  the  parallel  dielectric  component,  e||,  then  allows 
the  long  wavelength  whistler  wave  to  couple  to  the  short  wavelength  Bohm- 
Gross  mode  near  plasma  resonance.  The  whistler  wave  propagating  to  plasma 
resonance  then  mode  converts  into  a  short  wavelength  thermal  mode  that 
propagates  away  from  resonance  towards  decreasing  density. 

This  particular  mode  conversion  process  is  of  general  interest  because 
it  represents  a  channel  for  converting  wave  energy  into  fast  electrons. 

The  auroral  ionosphere  which  is  characterized  by  non- thermal  levels  of 
whistler  waves  and  field-aligned  electron  distributions  is  a  natural 
plasma  in  which  this  process  is  likely  to  occur.  In  the  auroral 
ionosphere  the  density  gradient  in  the  background  plasma  is  aligned  nearly 
along  the  magnetic  field  so  that  a  model  plasma  with  a  field  aligned 


density  gradient  is  appropriate  for  describing  this  environment.  Whistler 
mode  conversion  may  occur  naturally  in  the  auroral  ionosphere  or  be 
stimulated  artifically  by  launching  waves  from  polar  orbiting  spacecraft. 


The  main  purpose  of  this  study  is  to  determine  analytically  the 
detailed  structure  of  the  electric  fields  involved  in  the  mode  conversion 
process .  This  knowledge  can  then  be  used  in  later  computations  to 
determine  the  changes  in  the  electron  distribution  produced  by 
acceleration  from  the  mode  converted  whistler  wave. 

In  section  II  the  equation  governing  the  structure  of  the  electric 
potential  is  derived  and  an  integral  representation  of  the  general 
solutions  is  obtained.  Asymptotic  and  series  expressions  for  the  short 
wavelength  solutions  are  obtained  in  section  III  while  the  properties  of 
the  long  wavelength  solutions  are  examined  in  section  IV.  The  mode 
conversion  process  is  treated  in  section  V  and  the  structure  of  the  mode 
converted  electric  field  is  determined  in  section  VI.  Conclusions  are 
presented  in  section  VII. 

II .  Governing  Equation: 

The  equation  determining  the  structure  of  the  electric  potential  in 
the  plasma  is 

V-£-V  $  -  0,  (1) 

where  $  is  the  electric  potential  and  £  the  plasma  dielectric  tensor 
operator.  For  the  model  plasma  described  here,  the  electric  potential 
amplitude,  4> ,  varies  only  along  the  z-direction  and  $  can  be  written  in 
the  form 

$(x,z,t)  -  <^(z)exp[  i(k^x  -  wt)]  .  (2) 

Furthermore  the  plasma  dielectric  tensor  components  depend  only  upon  z 


and  in  the  cold  plasma  limit  are;  ~  ~  ^zy  ~  '^xx  “  ^yy  “ 

^xy  “  ■  i (^/‘*’)  ( £ 1 )  “  ■  ^yx •  ^'zz  *" 

A  plasma  with  a  nonzero  temperature  supports  modes  which  propagate 
both  along  (Bohm-Gross  modes)  and  perpendicular (Bernstein  modes)  to  the 
magnetic  field.  The  Bohm-Gross  modes  arise  from  thermal  corrections  to 
the  parallel  component  of  the  cold  plasma  dielectric,  £||  ,  while  the 
Bernstein  modes  arise  from  thermal  corrections  to  the  perpendicular 
component,  However,  since  we  are  considering  whistler  waves  which 

necessarily  have  w  <  0  we  need  not  consider  coupling  to  Bernstein  modes. 
In  addition,  we  consider  the  case  in  which  the  electron  gyroradius  (v^/fi) 
is  much  smaller  than  the  perpendicular  wave  length  so  that  thermal 
corrections  to  the  perpendicular  dielectric  component  are  negligible. 

The  lowest  order  thermal  corrections  to  the  parallel 
component  of  the  plasma  dielectric  cause  it  to  become  a  differential 
operator , 


where  v  is  the  thermal  velocity,  y(2T/m),  and  u>.p  is  the  electron 
plasma  frequency  which,  here,  is  a  function  of  z.  The  factor  7  has  the 
value  three  since  the  electron  motion  is  essentially  one  -  dimens ional .  The 
parallel  wave  number  usually  appearing  in  the  Bohm-Gross  dispersion 
relation  has  been  replaced  by  the  operator,  -i  d/dz ,  because  the  plasma 
is  nonuniform  in  the  z -direction. 

Since  the  mode  conversion  process  occurs  in  the  immediate  vicinity  of 
plasma  resonance  we  make  the  approximation  that  the  density  varies 
linearly  along  the  magnetic  field  with  scale  length,  L,  the  density  scale 
length  at  plasma  resonance.  The  plasma  frequency  can  then  be  written  as 


where  is  the  dimensionless  position  variable,  z/L.  Since  the  number 
density  in  the  plasma  is  always  positive,  must  be  restricted  to  the 
regime  less  than  unity  in  order  to  model  a  physical  plasma.  Using  (4) 
the  perpendicular  and  parallel  plasma  dielectric  components  become, 

c,  “  (?  -  Y2)/(1-Y2)  ;  c  =  j.  +  ^  £2^  ,  (5) 

where  Y  -  fi/w  and  -  yv^/(co^L^) .  Using  (2)  and  (5)  in  (1)  then 
results  in, 

-  p^<i>  -  0  ,  (6) 

which  determines  the  spatial  structure  of  the  electric  potential. 

Equation  (6)  has  the  same  form  as  the  equation  analyzed  by  Maggs ,  et.  al. 
(1984)  except  for  the  sign  of  the  term  proportional  to  <i> .  The  sign 
difference  arises  because  here  w  <  n,  while  the  earlier  work  treated  the 
case  w  >  Q.  In  (6)  the  parameters  and  ^2  depend  upon  the  perpendicular 
wave  number. 


Y2  -1 


'i  ^2^, 


(7) 


and  the  notation  denotes  the  derivative  with  respect  to  g- . 

The  parameters  and  ^2^  ere  real  and  positive  since  the 
electrostatic  mode  is  in  the  whistler  band  and  w  <  fl  so  that  Y  >  1.  The 
parameter  a  has  the  value  unity  for  a  plasma,  but  we  consider  general 
values  of  a  in  the  mathematical  analysis  of  Eq.  (6)  because  this  procedure 
is  necessary  to  obtain  the  structure  of  the  electric  field  in  the  vicinity 
of  plasma  resonance,  as  pointed  out  by  Rabenstein  (1958)  in  his  classic 
paper  on  mode  conversion.  The  parameter  -=  7v^/(w^L^)  is  the  square  of 


or  -  0) 


the  ratio  of  the  Debye  length  at  plasma  resonance  (i.e.  ,  u>  —  u> 
to  the  density  scale  length,  and  is  here  taken  to  be  a  small  parameter. 

If  is  set  equal  to  zero  a  second  order  differential  equation, 
referred  to  as  the  reduced  equation,  is  obtained  from  (6) .  The  reduced 
equation  determines  the  potential  structure  in  the  cold  plasma  limit. 
Transforming  the  independent  variable  of  the  reduced  equation  to  ^ 
and  substituting  the  function  -  u  exp(-^/2)  results  in 

Cu(2)  +  (a  -  e)u(l)  -(a/2  -  i^o)u  -  0  ,  (8) 

where  /Sq  “  Equation  (8)  has  the  form  of  Rummer's  differential 

equation,  and  linearly  independent  solutions  can  be  written  as  combinations 
of  the  Rummer  functions  M(a.,Q;;i^)  and  U(a.,a;i^),  where  a.  -  q/2  -  i^Q 
(Slater,  1970). 

In  addition  to  the  two  independent  solutions  to  the  reduced  equation, 
mother  pair  of  basic  solutions  exist  which  depend  upon  being  finite 
(i.e.,  f  0).  This  solution  pair  is  related  to  the  Bohm-Gross  modes  and 
thus  has  relatively  short  wave  length  in  comparison  with  the  solutions  of 
the  reduced  equation.  The  equation  describing  the  pure  thermal  modes  can 
be  obtained  from  (6)  by  setting  (  and  hence  equal  to  zero.  Doing 
this,  changing  the  independent  variable  to  7  -  and  writing  u(7) 

instead  of  4>ig-)  gives,  after  integration, 

u^^^  +  t7u(^)  +  (a  -  l)u  =  constant.  (9) 

For  the  case  a  -  1,  Eq.  (9)  becomes  the  inhomogeneous  (i.e.,  driven) 

Airy  equation  for  u^^)  ,  which  is  proportional  to  the  electric  field. 

The  short  wave  length  solutions  are  thus  related  to  Airy  functions 
and  the  integrals  of  Airy  functions. 

Since  the  coefficients  of  the  derivatives  in  (6)  are  all  constant  or 
linear  functions  of  the  independent  variable  ?,  general  solutions  can  be 


obtained  in  integral  form  by  using  the  Laplace  transform  technique 
(Coddington  and  Levinson,  1955).  The  solutions  then  have  the  form, 


S(?-) 


ds 


■s^ 


P(s) 


exp  (  -I  ®  Q(t)  dt]  , 


(10) 


P(t) 


where  P(s)  and  Q(s)  are  polynomials  obtained  from  (6)  by  replacing  i^(’^) 
by  (-s)^  to  obtain  an  expression  of  the  form,  P(s)?  +  Q(s) .  Using  this 
procedure  we  find  that  P(s)  -  s^  +  Q(s)  -  c^s*  -  as  -  The 

contours  denoted  by  C  in  (10)  are  chosen  so  that  the  bilinear  concomitant. 


e'^^  exp[  - r  Q(t)  dt)  , 
P(t) 


(11) 


which  is  essentially  the  boundary  term  obtained  after  an  integration  by 
parts,  is  zero  at  the  end  points  of  the  contour.  For  the  problem  under 
consideration  the  solutions  have  the  form 


(s  +  i/3  )“'*■(  s  -  i)S„)“'  exp[  -(c^s^/3  +  s?^)]  ds ,  (12) 


^  ^  , 


where , 

(13. a) 
(13. b) 

P  -  ;  Po  -  (13. c) 

The  contours  in  the  integral  representation  (12)  must  be  chosen  so  that 
the  condition 


a+  “  q/2  -  1  +  i)9  ;  a.  -•  q/2  -1  -  i^  , 

.20  3- 


(s  -  exp[  -(£2sV3  +  S?)]  -  0  ,  (14) 


is  satisfied  at  the  end  points. 


Since  and  a.  are  not  integers,  the  points  s  ■■  and  s  -  -1^2 
the  complex  s -plane  are  branch  points.  Therefore,  to  obtain  analytic 
solutions  using  the  integral  representation  it  is  necessary  to  cut  the 
integration  plane  with  branch  lines.  Since,  in  general,  f  a.  the  two 

branch  points  can  not  be  connected  by  a  single  branch  cut.  Thus  the 
integration  plane  is  cut  with  two  branch  lines  extending  from  the  branch 
points  to  infinity  parallel  to  the  real  axis  as  shown  in  Fig.  2.  In  the 
cut  plane  the  term  (s-i/92)*^’  is  single  valued  for  2mn  <  6  .<  2^  +  2m7r, 
where  8.  denotes  the  argument  of  s-i;92  and  m  is  integer.  The  term 
(s-t-i/92)^'*’  is  single  valued  for  2n5r  -  w  <  <  »r  +  2n»r,  where  denotes 

the  argument  of  s-t-i^2  ^  is  integer. 

Equation  (14)  can  be  satisfied  by  contours  that  end  at  the  points 
s  -  or  s  -  •i-^2'  Pi^ovided  that  Re(Q/2  ±  i/9)  >0,  or  by  contours  that 
proceed  to  infinity  along  rays  such  that  exp(-£^s®/3)  vanishes.  It  is 
convenient  to  define  two  basic  classes  of  solutions  with  members  denoted 
by  and  .  The  solutions  Aj  are  obtained  from  the  integral 
representation  (12)  by  using  contours  that  begin  and  end  at  infinity.  The 
solutions  Bj  are  obtained  by  using  contours  with  at  least  one  end  point  at 
s  -  i/32  •  Some  of  these  contours  are  illustrated  in  Fie.  2.  The 

solutions  B2,  B^,  and  Bq  are  closely  related  to  the  solutions  of  the 
reduced  equation  given  in  (8).  Asjnnptotically  B2  and  B^  represent  long 
wavelength  propagating  waves.  The  solution  A^^  is  related  to  the  /  0 
class  of  solutions  which  satisfy  Eq.  (9).  Asymptotically  A^  represents  a 
short  wavelength  propagating  wave. 


III.  ohort  Wave  Length  Solutions: 


In  this  section  we  investigate  the  behavior  of  the  solutions  that 

asymptotically  correspond  to  the  short  wavelength  solutions  of  Eq .  (9). 

We  evaluate  the  leading  term  of  the  asymptotic  series  for  the  solutions 

which  is  useful  for  determining  their  behavior  at  large  .  We  then  obtain 

a  series  representation  for  A.  which  can  be  used  to  evaluate  these 

J 

solutions  for  small  'i- .  The  asymptotic  expression  is  needed  to  identify 
the  thermal  mode  produced  in  the  mode  conversion  process  while  the  series 
expansion  is  needed  to  evaluate  the  electric  field  structure  near  the 
origin  (i.e.,  plasma  resonance). 

The  contours  describing  the  solutions  A^  have  both  end  points  at 
infinity.  An  example  of  the  contour  for  the  solution  labeled  A^  is  shown 
in  Fig.  2.  In  order  for  the  condition  (14)  to  be  satisfied  at  large 
Isl,  the  term  c^s^/3  must  have  a  positive  real  part,  i.e, 

27r(2-j)/3  -7r/6  <  arg(s)  <  »r/6  +  2;r(2-J)/3  .  (15) 

The  integer  j  in  (15)  labels  the  open  sector  in  the  s-plane  of  angular 
width  ff/3  centered  about  the  angle  27r/3  for  j-1,  0  for  j-2  ,  and  -27r/3  for 
j-3.  The  contour  for  the  function  A^  has  its  end  points  at  infinity  in 
the  two  sectors  other  than  sector  j.  For  example,  the  contour  for  the 
function  A^  proceeds  to  infinity  in  the  sectors  labeled  by  j  -  2  and  3, 

In  order  to  study  the  properties  of  the  functions  A^  it  is 
convenient  to  make  the  transformation  of  variables 

a  -  ^  _  £‘2/3^  ;  Op  -  i£2/3^^  (2^6) 

With  this  transformation,  the  integral  representation  (12)  has  the  form 


(ct-CT(,)“' (a+Op)®'''  exp[  -((7^/3  +  077)]  da  ,  (17) 

C(Aj) 


A  j  ( 17 )  -  e 


-2(a-l)/3 


where  C(Aj)  is  the  image  of  the  contour  for  under  the  transformation 
(16)  . 

We  first  investigate  the  asymptotic  properties  of  the  functions 

by  assuming  that  is  small  and  the  magnitude  of  ^  is  large, 

specifically,  1^1  »  ■  The  saddle  point  method  of  integration  can  be 

used  to  evaluate  (17)  under  these  asumptions.  Defining  f(o)  -  ct^/3  +  ar/ , 

the  saddle  points  are  determined  by  the  condition  f (a)  -  0.  Thus  the 

integrand  of  (17)  has  two  saddle  points  at  The  condition  of 

large  I  I  insures  that  cr^  »  ,  i.e.,  the  saddle  points  are  far  removed 

from  the  singularities  in  the  integrand  of  (17).  Demanding  that  a  saddle 

point  lie  between  the  Stokes  lines  (which  are  located  at  arg(s)  -  0,  2n/3, 

and  47r/3)  places  the  following  restrictions  on  the  argument  ranges  of  the 

functions  A.  , 

J 

2ir<2-j)/3  +  2k/3  <  arg(a^)  <  ^n/3  +  2n(2-j)/3  .  (18) 

For  the  saddle  point  a  -  (18)  gives  the  argument  ranges 

47r(2-j)/3  +  7r/3  <  0  <  5n/3  +  47r(2-j)/3  ,  (19. a) 

while  the  saddle  point  a  -  gives  the  argument  ranges 

Un(2-j)/3  -  5k/3  <  8  <  -k/3  +  47r(2-j)/3  ,  (19. b) 

where  6  -  arg(?-)  .  The  argument  ranges  for  the  two  values  of  the  saddle 
point  are  illustrated  in  Fig.  3. 

The  direction  of  the  path  of  steepest  descents  through  the  saddle 
point  is  obtained  by  assuming  f((7)  has  the  form,  -  f "  (o^)  (a-o^)  2/2  , 
where  is  real  and  positive  (i.e.,  has  argument  zero).  The  sign  of  the 
leading  term  in  the  saddle  point  method  of  integration  is  determined  by 
choosing  the  sign  of  x  -  ±(o'-<7q )  (f "  (og)  ^  ^ 2 )  such  that  x  is  positive  along 
the  direction  of  the  contour  of  integration.  Thus,  if  the  contour  of 


integration  passes  through  the  saddle  point  in  the  same  (opposite) 
direction  as  the  path  of  steepest  descents  the  plus  (minus)  sign  is 


chosen.  Since  f"(ag)/2  -  ,  the  positive  direction  of  the  path  of 

steepest  descents  has  arg(a-a^)  -  -arg [( ir? ^  ^ ^ /2  for  the  saddle  point 

and  arg(a-a^)  -  -arg  [( -ir?^ ''^ )  i  /  2  j /2  for  <7^  -  .  The  leading 

term  in  the  asymptotic  expression  for  then  gives,  for  the  saddle  point 

ag  = 


Aj  (ry)  “  (-)^  Pj  y^r  (irj^  ^  ^  £'2(q-1)/3  g  ^ 


(20. a) 


where  P2  -  p^  “  1.  and  p^  -  .  The  factor  p^  in  (20)  arises 

because  the  saddle  point  for  the  solution  lies  in  the  Riemann  sheet 
where  the  argument  of  the  term  (a-(7Q)®+  is  In  larger  than  in  the 
principal  sheet.  For  the  saddle  point 

(20. b) 

In  (20. a)  and  (20. b)  the  variable  f  is  defined  as  f  •=  2rj®'^/3  . 


Aj  (»j)  —  (-)'^'''^  Pj  Jn  ( -  if;  ^  £'2(a-l)/3  g^f_ 


The  asymptotic  behavior  of  Aj  as  given  by  (20. a)  and  (20. b)  contains 
the  exponential  term,  exp(±if),  where  f  When  C  is  imaginary, 

the  argument  of  the  exponential  term  (i.e.,  ±iC)  is  purely  real  and  Aj  is 
either  spatially  growing  or  decaying.  For  example,  with  arg(?-)“jr/3 , 

C  -  il  and  ,  as  illustrated  in  Fig.  3,  A2 ,  which  is  asymptotically 

proportional  to  exp(-iC)  as  given  by  (20. a),  represents  a  spatially  growing 
function  while  (20. b)  indicates  A^  represents  a  spatially  decaying 
function. 

When  the  variable  ^  is  real  the  exponential  term  in  (20.a,b)  is 
oscillatory  and  Aj  represents  a  propagating  wave.  The  assumed  temporal 
dependence  exp(-iwt)  determines  whether  the  phase  velocity  of  the  wave  is 
positive  or  negative.  For  example,  with  arg  g  -  0,  C  is  real  and  positive 
and,  from  (19. a)  and  (20. a),  we  see  that  A^  is  defined  for  arg  g  -  0  and 


varies  as  exp(-iO.  while  (19. b)  and  (20. b)  indicate  is  defined  for  this 
argument  value  but  varies  as  exp(iO-  The  phase  of  then  varies  as 

(f-cjt)  so  that  A^  asymptotically  represents  a  positive  phase  velocity  wave 
(i.e.,  a  wave  propagating  left  to  right)  for  arg  ?•  -  0.  Similarly  A^  , 
for  arg  ■=  0,  represents  a  right  to  left  propagating,  or  negative  phase 
velocity,  wave.  The  asymptotic  behavior  of  A^  at  various  values  of  the 
argument  of  ^  is  indicated  in  Fig.  3. 

The  solutions  A^  can  be  evaluated  for  small  values  of  the  independent 
variable  ^  (i.e.,  near  plasma  resonance)  by  deriving  a  series  expression 
from  the  integral  representation  (17) .  This  can  be  accomplished  by 
expanding  the  term  (a+aQ)°‘'^  in  the  integrand  of  (17)  in  a  series 

in  Oq/o  and  integrating  term  by  term.  Carrying  out  this  procedure  (the 
details  of  the  calculation  are  given  in  Appendix  I)  gives  the  following 


series  representation  of  Aj  ,  valid  in  the  finite  -plane. 


Aj(,)  -  . 


•2(a-l)/3 


Pj  L  L  ^m  Ij(c*-m.k)  r?  , 

km 


where 

Ij(a-m,k)  -  27ri  exp  (,ri(l-p)  (7-2j )]  /  [3^  r(p)]  ,  (22) 

and  p  -  (4-a-m-k)/3.  The  phase  factors  pj  are  the  same  as  given 
previously  for  the  asymptotic  expressions  (20. a, b).  The  restricted 
argument  ranges  given  in  (19. a, b)  that  apply  to  the  asymptotic  expressions 
(20. a, b)  do  not  apply  to  the  series  expression  (21)  which  is  valid  for  all 
values  of  arg(r?)  .  Since  Uq  -  assumed,  for  a  plasma  with 


a  gentle  density  gradient,  to  be  small  because  e*  <  1.  With  Uq  small,  A^ 
can  be  approximated  from  (21)  by  using  only  the  m  —  0  term. 


'  liF.  nr:  ir  /  WV  W\J  wv 


IV .  Long  Wave  Length  Solutions: 

In  this  section  we  obtain  approximate  expressions  for  the  solutions  B. 
and  Bq  in  terms  of  the  well  known  Rummer  functions  M  and  U  (Slater,  1970). 
The  leading  term  of  these  approximations  is  related  to  the  cold  plasma 
electrostatic  whistler.  The  expressions  obtained  here  are  used  to 
identify  the  solutions  that  as)rmptotically  represent  whistler  waves  in  the 
mode  conversion  solution  given  in  Sec.  V. 

The  contours  for  the  solutions  Bj  and  B^  have  at  least  one  end  point 
at  the  branch  point  located  at  s  -  obtain  an  expression  for  B^ 

and  Bq  in  terms  of  familiar  functions  it  is  useful  to  express  the 
integrand  of  (12)  as  a  Taylor  series  in 

Bj(  £2)  _  B^(0)  +  £2  3^2Bj(0)  +  ...  ,  (23) 

where  B^ (0)  is  obtained  by  setting  £2  -  0  in  (12), 


i,<0)  -  (s  +  <3S  ,  (24) 


where  a_|_  is  evaluated  at  £2  —  0  ,  namely 

a+  -  a/2  ±  i^o  -  1  .  (25) 

In  (23)  and  (24)  j  can  have  the  values  0,  1,  2,  or  3.  Since  Bj (0)  is 
obtained  by  setting  £2  -  0  it  is  clearly  a  solution  of  the  reduced 
equation  (i.e.,  (6)  with  £2  -  0).  The  c2  expansion  given  in  (23)  thus 
represents  successive  thermal  corrections  to  the  cold  plasma  whistler 
mode . 

The  solution  B^  is  obtained  from  the  contour  which  starts  at  s  —  -^^2’ 
circles  the  point  s  -  in  the  counter-clockwise  direction  and  returns 
to  the  point  s  -  -i^2>  shown  in  Fig.  2.  Changing  the  variable  of 


integration  in  (24)  to  s  -  gives 

Bq(0)  =  (i/92)  J  (P-1)  ■  (P+1)  exp ( - ip)92?- )  •  (26) 

In  (26)  the  limits  of  integration  with  the  notation  (1+)  denotes  that  the 
contour  starts  at  p  -  -1,  circles  the  point  p  -  1  in  the  counter-clockwise 
direction,  and  returns  to  the  starting  point  p  -=  -1.  Since  the  integral 
representation  for  the  Rummer  function  M(b,c;x)  for  Re(b)  >  0  is  (Erdelyi, 
1953), 

r(f'^)  vt  b“l  c-b-l 

M(b,c:x)  -  r(c)r(b-c-H)  e^’'  t""  ^  (t-l)"^  dt  ,  (27) 

2jrir(b)  -^0 

Bq  can  be  obtained  in  terms  of  M(a,b;x)  by  changing  the  variable  of 
integration  in  (26)  to  t  -  (p-Hl)/2,  giving 

Bq(0)  -  (2i^2>“'^  exp(-i2/92^t)  dt.  (28) 

Choosing  b  -  ^^0  “  a  ■  a  and  c  -  a  in  (27)  then  gives, 

B  (0)  -  (2i^2)“'^  27ri  r(Q-a)  K{a-a,a; -2ip  g.) .  (29) 

r(o)r(l-a) 

where  a  -  a/2  -  i/Sg  -  a.  +1.  The  expression  (29)  is  valid  for  all 
values  of  arg(?-)  .  The  higher  order  corrections  to  Bq  in  the  expansion 
(23)  can  be  shown  to  depend  upon  M  and  its  derivatives  with  respect  to 
(Maggs,  et.al.,  1984),  but  we  do  not  attempt  this  here. 

The  contour  for  the  solution  (0)  extends  from  s  -  -^^2  infinity 
in  the  sector  labeled  'j'  in  Figure  1,  passing  above  the  point  s  —  i;32  • 

At  large  I  si  the  integrand  in  (24)  varies  as  s“*2  .  In  order  for  the 


integrand  to  be  bounded  along  the  contour  for  q  >  1,  it  is  necessary  to 


require  that  ^  be  nonzero,  i.e.  I?- 1  >  l^pl  >  0  where  ?-q  is  arbitrary  but 


small  and  non  zero.  Thus  for  the  solutions  (j  -  1,2,3)  the  expansion 
(23)  is  valid  only  for  nonzero  values  of  ?• .  This  feature  is  evident  in  the 
Kummer  function  U(b,c;x)  which  is  singular  at  the  origin  for  Re(c)  >  1. 

In  addition  to  the  requirement  I?- 1  >  Ij-qI,  the  integrand  in  (24)  is 
bounded  only  when  Re(s?-)  >  0,  i.e., 

-jr/2  <  arg(s?^)  <  ir/2  .  (30) 

The  requirement  (30)  combined  with  (15)  gives  the  argument  ranges  over 
which  the  functions  (0)  are  defined  by  the  integral  representation  (24) , 
2ir(j-2)/3  -  2V3  <  arg(?^)  <  2jr/3  +  2»r(j-2)/3  .  (31) 

The  argument  ranges  given  by  (31)  are  compared  to  those  for  the  functions 
Aj  as  defined  by  (19. a, b)  in  Figure  4.  Using  the  transformation 


s  -  2i)32(t+l/2)  , 


in  (24)  results  in 


B.(0)  -  (2i;92)“'^  exp(-2i^2^t)  dt  ,  (33) 


where  the  contour  in  (33)  extends  to  infinity  along  a  ray  with 
arg(t)  -  B. 

The  integral  representation  of  the  Kummer  function  U(b,c;x) 


is  (Erdelyi,  1953) 

U(b,c;x)  -  r(l-b)  e 
2»ri 


-i^b  r(o+)  ^b-i  (i.,t)C-b-i 


exp(-xt)  dt.  (34) 


In  (34)  the  contour  starts  at  infinity  with  arg(t)  -  B,  circles  the 
origin  in  the  counter-clockwise  direction  and  returns  to  infinity  with 
arg(t)  -  ^  +  27r.  The  integral  in  (34)  can  be  divided  into  three  separate 
integrals , 


WTfWH  ^ 


mr-  -K" 


r(0+)  r-1  r(0+)  p»e 

ooe^^  ooe^^  -1  -1 


i0+2n 


.  I2  -H  I3  . 


The  integral  on  the  right  hand  side  of  (35)  is  like  except  for  the 


direction  along  the  contour  and  the  argument  of  the  variable  of 


integration.  The  argument  of  the  variable  of  integration  in  is  2n 


larger  than  in  so  that  the  term  in  is  exp(i2»rp)times  the  same 

term  in  .  Since  the  argument  difference  of  2n  does  not  alter  the  other 


terms  in  the  integrand  of  (34),  is  simply  -exp(i27rp)  times  .  The 


integral  in  (35)  can  be  related  to  e  ^2^  M(Q-a,Q; -2i/92?-)  by  changing 


the  variable  of  integration  from  t  to  t  -  t+1.  Applying  (35)  to  (33)  and 


using  (34)  then  gives 


B  (0)  -  (2i/3  )^-^  TT  e~^^2^  (  U(a,Q;2i^  ^) 

sin(7ra)r(l-a) 


e  r(a-a)M(a.a:2i5^g^); 
r(a)  ^ 

(36) 


where  a  -  a/2  -  i/Sg .  Even  though  the  subscript  j  does  not  appear 


explicitly  on  the  right  hand  side  of  (36)  it  is  present  implicitly 


because  of  the  argument  range  restrictions  given  by  (31).  Since 


U(a,a;2i^2^)  ^  multivalued  function  the  asymptotic  form  for  ,  as 


obtained  from  (36),  depends  upon  the  value  of  arg(?)  and  thus  on  the 


value  of  j . 


Equation  (36)  can  be  used  to  obtain  asymptotic  expressions  for  Bj 


that  can  be  related  directly  to  whistler  waves  in  the  cold  plasma  limit 


when  the  physical  value  unity  is  used  for  the  parameter  a  (i.e.  a-1) . 


For  example,  defining  w  -  2i^^^  and  taking  arg(?)-  0  gives  arg(w)  -  ff/2 


and  for  this  argument  value  the  asymptotic  values  of  U(a,l;w)  and 


M(a,l:w)  are 


U(a,l;w)  ~  w‘®  , 


(37. a) 


M(a,l;w)  -  w'^  +  ^a-l  gW 

r(l-a)  r(a) 


(37. b) 


From  (31)  and  Figure  4  we  see  that  the  argument  range  of  B2  contains 


arg  ?•  -  0  so  that  using  (37. a, b)  in  (36)  gives 


6^(0)  ~  -r(l-a)  e-i-^a  „a-l  eW/2 


With  the  chosen  time  dependence  exp(-ia)t),  the  phase  of  ,  as  given  by 


(38)  ,  varies  as  2/32?--cJt  so  that  B2  is  a  left  to  right  traveling  wave  for 


arg  ?•  -  0,  because  for  w  <  O  (i.e.,  electrostatic  whistlers)  the  phase 


velocity  and  group  velocity  are  in  the  same  direction  (forward  wave). 


As  a  further  example,  consider  the  asymptotic  behavior  of  B^  at 


arg  ?■  -  TT,  i.e.  arg  w  -  3jr/2.  The  asymptotic  values  of  U(a,l;w)  and 


M(a,l;w)  at  arg(w)  -  3jr/2  are, 


U(a,l:w)  -  w'^  .  27ri  w*^  .  2ni  w^'^  e^  , 

r(a)r(l-a)  r(a)r(a) 


(39. a) 


M(a,l;w)  -  e'^^^  w~^  4.  w^'^  e^  . 

r(l-a)  r(a) 


(39. b) 


where  the  principal  branch  of  the  multivalued  function  U(a,l;w)  is  taken 


as  -n  <  arg(w)  <  jr.  From  (31)  we  recognize  that  the  argument  range  of  B^ 


contains  the  value  arg  ^  ~  re  so  that  using  (39. a, b)  in  (36)  gives 


B2(0)  ~  -r(l-a)  w^-l  e^’^a  eW2  . 


Since  the  argument  of  ^  is  jr,  the  phase  of  B^  varies  as  For 


constant  phase  then,  \^\  decreases  as  t  increases,  and  since  ^  is  negative 


5^  represents  a  left  to  right  propagating  wave  for  arg(?)—  tt. 


direction  of  propagation  for  Bj  with  a  -  1  is  shown  for  selected  values 


of  arg  ^  in  Fig.  4. 


The  asymptotic  expressions  (38)  and  (40)  are  consistent  with  the  WKB 


TW 


form  of  the  whistler  in  the  cold  plasma  limit.  Using  Eq .  (6)  with  -  0 


and  the  fact  that  the  WKB  solutions  vary  as  exp  (ij  k||d?^)  ,  gives  the  WKB 


dispersion  relation, 


3-  k||2  +  iQk||  -  0  , 


where  the  derivative  of  k  has  been  assumed  negligible.  From  (41)  we 


obtain 


k.,  =  i(a/2  +  i-fio)/?-  +  /3  , 


where  terms  of  order  1/3-^  have  been  ignored  since  3^  is  assumed  large. 


Using  (42)  then  gives 


iJ  k  d#  -  -(a-a)ln^  +  3  , 


so  that  the  WKB  solutions  are  proportional  to  (3)^'°‘ex-p(.ip^3) .  For 


a  =  1  this  behavior  is  the  same  as  obtained  in  ('3fi)  and  (40)  for  B„ 


and  B^ . 


The  functions  (0)  (with  j  -  1,2,3)  are  singular  at  the  origin 


because  of  the  singularity  in  the  function  U(a,a:  21^2?’)  •  a  -  1  the 


function  U  diverges  as  In^  as  3  approaches  zero.  Thus  near  the  origin  the 


expansion  of  (23)  is  not  useful  for  evaluating  B^  However,  the 


behavior  of  the  solutions  B^  near  the  origin  can  be  determined  by 


expressing  them  in  terms  of  the  solutions  B^  and  as  discussed  in 


section  VI . 


V .  Mode  Conversion: 

In  this  section  we  consider  the  case  of  an  electrostatic  wave 
propagating  in  the  whistler  band  incident  upon  plasma  resonance  from 
inside  the  plasma.  In  the  whistler  band  the  phase  velocity  is  in  the  same 
direction  as  the  group  velocity  (i.e.,  the  whistler  is  a  forward  wave)  so 
that,  in  the  model  plasma,  the  incident  whistler  is  a  positive  phase 
velocity  wave  approaching  the  origin  along  the  negative  real  ?--axis.  Thus 
we  look  for  a  solution  at  negative  values  of  ? ,  i.e.  arg(?)  -  jr  which 
contains  a  part  corresponding  to  a  long  wavelength  incident  wave  (i.e., 
propagating  from  left  to  right)  together,  perhaps,  with  a  reflected  wave. 
On  the  other  side  of  plasma  resonance  (i.e.,  arg(?^)-0)  the  solution  must 
be  bounded  at  infinity  and  be  consistent  with  energy  transport  away  from 
the  origin.  Thus  for  arg  g-  -  0  the  solution  must  contain  left  to  right 
propagating  waves  only. 

From  the  discussion  in  section  IV,  a  candidate  solution  for  the 
incident  wave  is  since  it  is  long  wavelength  and  propagates  left  to 
right  for  arg  -  n.  From  Fig.  4  we  see,  however,  that  the  argument 
range  of  does  not  contain  arg  ^  -  0  so  that  we  can  not  directly 
determine  the  properties  of  B^  on  the  other  side  of  plasma  resonance.  To 
evaluate  B^  for  arg  ?•  -  0  we  must  analytically  continue  B^  into  this 
argument  regime.  Refering  to  Fig.  2  we  note  that  the  cut  integration 
plane  is  free  from  singularities  in  the  regions  between  the  contours 
defining  the  solutions  B^ ,  B2  and  .  The  solution  B^  can  then  be 
related  to  the  solutions  B2  and  A^^  by  extending  the  contour  for  B^  to 
infinity  in  sector  2  and  applying  Cauchy's  Theorem.  This  procedure  yields 


the  expression 


B,  -  B„  -  e-i2ira.  a,  . 


whi^i'  analytically  continues  into  the  desired  argument  range  since  as 
shown  in  Fig.  4  the  argument  ranges  of  and  both  contain  arg(?^)  -  0. 

The  factor  e'^^Tra.  (44)  occurs  because  part  of  the  distorted  contour 
for  B^  (the  dashed  portion  in  Fig.  2)  lies  on  the  Riemann  sheet  for  which 
the  phase  of  the  term  (a-Og)®-  in  (17)  is  2»r  less  than  it  is  on  the 
principal  sheet  in  which  the  defining  contour  for  lies.  The  fact  that 
the  analytic  continuation  involves  the  short  wave  length  mode  represented 
by  A^  means  that  mode  conversion  must  occur.  Thus  the  mathematical 
procedure  of  analytic  continuation  corresponds  to  the  physical  process  of 
mode  conversion. 

B^  analytically  continues  into  a  physically  acceptable  solution 
because,  for  arg(^)-0,  (38)  indicates  that  asymptotically  (I?|-»  ®,  arg(?-)-0) 
B2  corresponds  to  a  left  to  right  propagating  wave,  while,  from  Fig.  3,  we 
note  that  A^  also  corresponds  to  a  left  to  right  propagating  wave  at  arg  i- 
-  0.  Thus  both  waves  transport  energy  away  from  plasma  resonance.  We 
find  then  that  a  left-to-right  propagating,  purely  long  wavelength  mode 
(represented  by  B^)  incident  upon  plasma  reson;.nce  from  inside  the  plasma 
gives  rise  to  both  a  long  wave  length  transmitted  wave  (represented  by  B2) 
and  a  short  wave  length  transmitted  wave  (represented  by  A^^) .  The 
transmission  coefficient  for  the  cold  wave  is  or  ie*’^^o.  The 

reflection  coefficient  is  zero,  i.e. ,  there  is  no  reflected  wave.  These 
results  are  the  same  as  obtained  by  Bafios,  et.  al.  (1986),  who  considered 
only  the  purely  cold  case  represented  by  the  reduced  equation  (i.e.,  (6) 
with  -  0) . 

In  obtaining  the  asymptotic  expressions  for  B^ ,  B2 ,  and  A^  we  have 
assumed  that  the  magnitude  of  the  function  arguments  is  large.  While  this 
assumption  is  consistent  with  the  model  plasma  for  negative  values  of  , 
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it  is  not  for  positive  values  of  ^  when  considering  the  function  . 

This  is  because  the  region  ?  >  1  corresponds  to  a  fictitious  negative 
value  of  plasma  density  arising  from  the  approximation,  1  -  -  z/L, 

in  the  model.  Thus  while  the  argument  of  the  function  f 
can  be  assumed  large  for  values  of  ?  <  1  because  <  1,  the  argument  of 
the  function  B2  can  not.  However,  an  appropriate  expression  for  the 
function  B2  can  be  obtained  by  assuming  that  the  magnitude  of  the 
coefficient  a  in  the  Kummer  functions  U(a,a:w)  and  M(a,a;w)  is  large  in 
comparison  to  a  and  the  magnitude  of  the  indep'^nent  variable  w.  In  terms  of 
physical  parameters  this  requirement  becomes, 

1^1  <  .  (^5) 

Since,  for  waves  propagating  in  the  whistler  band,  co  <  n  the  requirement 
(45)  can  be  satisfied  for  values  of  p  with  <  ?•  <  1  ,  where  again  is 
arbitrary  but  small  and  non  zero.  Then,  for  arg(^)  -  0  and  large  lal, 
U(a,l;w)  and  M(a,l:w)  have  the  as3mptotic  behavior  (Slater,  1970) 
e''^/^  U(a,l;w)  ~  r(l-a)  [cos(wa)Jo  (ix)  •  sin(;ra)Yo  (ix)  ]  ,  (46. a) 

e-V2  M(a,l:w)  -  Jodx)  .  (46. b) 

where  Jq  and  Yg  are  zero  order  Bessel  functions  of  the  first  and  second 
kind,  respectively,  and  x  “  2(2/3o/32^)  ’ -  Using  (46. a, b)  in  (36)  we  then 
obtain 

B2(0)  -  i7r[Jo(ix)  +  iYo(ix)]  -  2Ko(x)  ,  (47) 

where  Kg  is  the  modified  Bessel  function.  From  (47)  we  see  that,  for 
arg(?)  -  0,  B2  represents  a  decaying  function  for  increasing  1^1. 

Thus,  a  long  wavelength  whistler  incident  upon  plasma  resonance  from 
inside  the  plasma  gives  rise,  beyond  plasma  resonance,  to  a  short  wavelength 
Bohm-Gross  mode  propagating  towards  decreasing  density  and  a  decaying  long 
wavelength  structure.  No  energy  is  reflected  back  into  the  plasma  from 


w 
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VI .  Electric  Field  Structure  Near  Plasma  Resonance 

Having  established  that  a  whistler  wave  incident  upon  plasma 
resonance  converts,  without  reflection,  into  a  short  wavelength  thermal 
mode,  we  now  determine  the  structure  of  the  electric  field  near  plasma 
resonance  arising  from  the  mode  conversion  process.  We  then  relate  the 
amplitude  of  the  field  at  the  origin  to  the  asymptotic  amplitude  of  the 
incident  whistler  wave. 

The  structure  of  the  electric  field  of  the  incident  cold  mode  in 
the  vicinity  of  plasma  resonance  can  not  be  obtained  from  expressions 
such  as  (23)  because  the  function  singular  at  the  origin.  The 

structure  can  be  found,  however,  by  relating  the  solution  to  the 
solutions  Aj  and  Bq  which  are  not  singular  at  the  origin.  In  addition, 
this  relation  gives  an  expression  for  B^  valid  for  all  values  of  arg(?^) 
when  the  series  expressions  for  the  functions  are  used. 

The  relation  useful  for  evaluating  the  behavior  of  B^  near  the  origin 
is  established  by  deforming  the  contours  which  define  the  solutions  and 
A2 .  As  shown  in  Fig.  5.  a,  B^  can  be  related  to  B^^ ,  A^  and  A^^  by  extending 
the  contour  defining  B^  to  infinity  in  sectors  1  and  2.  This  procedure  is 
permitted  because  the  integrand  in  (12) ,  the  general  expression  for  the 
solutions,  is  analytic  in  the  cut  s-plane  over  the  region  of  distortion. 
Then  by  Cauchy's  theorem  we  obtain  the  relation 

Bj  -  B^  -  Aj  -  e-i2’ra.  (48) 

The  factor  e'^^jra.  multiplying  Aj^  arises  because  the  argument  of  the  term 
(s  -  in  (12)  is  2n  less  than  over  the  basic  defining  contour  for  A^ 

shown  in  Fig.  2.  Likewise,  the  contour  defining  can  be  pulled  into  the 
point  s  -  ‘i^2  coincide  with  the  contours  for  B^  and  the  combination 

gi2ffa+  [gi27ra.  shown  in  Fig.  5.b,  giving 


■ 
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-  -B^  -t-  ei2’ra+  [ei27ra.  5^  +  (49) 

The  relations  (48)  and  (49)  can  then  be  combined  to  give 

B3(?)  -  [ei^iTQ  .  1]-1[  A2  +  e-i^Tra.  -  ei-2wa+  j  (50) 

The  expression  (50)  for  B^(^)  is  valid  for  all  arg(?-)  when  the  series 
expressions  (21)  for  A^  are  used,  so  that  (50)  represents  the  analytic 
continuation  of  B^  into  all  arg(?)  values  for  the  finite  ?^-plane. 

Note  from  (50)  that  the  term  outside  the  braces  diverges  as  q 
approaches  integer  values,  so  that  it  must  be  the  case 

Lim  [  A^  +  e'i^Tra.  a^  A^  -  ei^’ra^.  ]  -  0  .  (51) 

Q  -»  1 

Using  L'Hopital's  Rule  to  evaluate  (50)  we  then  obtain 

B3  -  +  Oa(e-i2^a.A^)  +  a^A^  -  ac»(ei27ra+BQ)  j  | J  .  (52) 

2?ri 

Note  that  to  obtain  (52)  it  is  necessary  to  treat  a,  the  coefficient  of  the 
first  derivative  term  in  (6),  as  a  general  parameter  as  was  first  pointed 
out  by  Rabenstein  (1958)  in  his  extension  of  the  work  by  Wasow  (1953). 

The  function  proportional  to  the  electric  field,  B^'  (where  ' 
denotes  derivative  with  respect  to  argument) ,  is  given  by 

B3'  -  _1_  (  d^A^’  +  a„(e-i2»ra.A^')  -h  aoA3')  ,  (53) 

27ri 

where  the  term  proportional  to  the  derivative  of  has  been  ignored  in 

2  /s 

(53)  because  it  is  smaller,  by  a  factor  of  e  ,  than  the  terms 
proportional  to  the  derivatives  of  the  functions  A^  .  Assuming  that 
(7g  -  <  1,  the  m  -  0  term  of  the  series  expansion  (21)  can  be  used 

to  evaluate  (53)  to  obtain  B3'  in  terms  of  familiar  functions, 

B3'  -  n  [  Gi(-rr)  +  iAi{-v)  ]  ■  (54) 

The  details  of  this  calculation  are  given  in  Appendix  II.  In  (54),  Gi(-p) 


denotes  a  solution  of  the  inhomogeneous,  or  driven,  Airy  equation 
and  Ai(-r;)  is  a  solution  to  the  homogeneous  Airy  equation  (Antosiewicz , 


1970) .  The  asymptotic  form  of  the  functions  in  (54)  are 
Ai(-f?)  -  sin(C  +  7r/4)  , 

and  Gi(-rj)  -  cos(C  +  n/U)  - 


(55. a) 


(55. b) 


where  C  -  2t7®''^/3.  From  (55. a, b)  it  is  clear  that  the  electric  field 
structure  of  the  mode  converted  wave  consists  of  a  short  scale  length 
wave  propagating  down  the  density  gradient  away  from  plasma  resonance 
together  with  a  field  that  falls  off  slowly  as  1/r?. 

The  electric  field  amplitude  at  the  origin  can  be  related  to  the 
asymptotic  whistler  amplitude  by  obtaining  an  approximation  to  Eq.  (6) 
valid  near  the  origin  for  0=1.  With  a  -  1,  (6)  can  be  integrated  once 


to  obtain 


ds-  +  C  , 


where  C  is  a  constant  of  integration.  Now  consider  the  mode  conversion 
solution  which  satisfies  (56).  Near  the  origin  8^(0),  which  is  the 
leading  term  in  the  expansion  of  is  singular  for  small  9-,  diverging 
as  In(^).  For  the  purpose  of  estimating  the  relative  magnitudes  of  the 
terms  in  (56)  for  small  ?• ,  B^  (and,  thus  4>)  can  be  approximated  by  ln(?-)  . 

Since  is  proportional  to  is  of  order  unity,  while 

£^(^(3)  is  proportional  to  The  terms  on  the  left  hand  side  of  (56) 

are  then  of  order  unity,  or  larger,  for  small  ?• .  On  the  other  hand,  the 
integral  term  on  the  right  hand  side  of  (56)  approaches  zero.  Evaluating 
this  term  we  find, 

f  -  ^2^  if^)ln(?^)d?  -  .  (57) 

Since  and  are  continuous  and  have  the  limit  zero  as  ^  goes  to 

zero,  the  integral  term  on  the  right  hand  side  of  (56)  can  be  made  as 


small  as  desired  by  chosing  ?•  small  enough.  Thus  for  ?•  <  where  is 


such  that 

?o(ln(?o)-l)  -  «  Min(C,l)  ,  (58) 

(56)  can  be  approximated  by 

£2  ^(3)  +  ^^(1)  -  c  .  (59) 

By  setting  (/i(l)  -  E,  the  z-component  of  the  electric  field,  and  using 

z  -  gL,  (59)  can  be  written  in  the  form 

7  _dl  E  +  z  E  -  Eq  .  (60) 

“^pe^  dz^  L 

Equation  (60)  describes  plasma  resonance  driven  by  a  constant  electric 
field  of  amplitude  Eq .  Thus  we  interpret  the  constant  C  in  (59)  as  an 
effective  external  electric  field  strength  at  the  origin,  Eg,  which 
excites  the  Bohm-Gross  mode.  The  solution  to  (59)  is  then 

^(1)  -  -JT  Eo  £-2/3  [Gi(-.7)  +  iAi(-r?)]  .  (61) 

Interpreting  (^(^)  as  an  approximation  to  the  solution  B^',  i.e., 

,^(1)  c  S^B^  where  C  is  a  constant  of  proportionality,  we  obtain 

C  Sj,  B3  -  -TT  £-2/3  Eo  lGi(-»7)  +  iAi(-r;)]  .  (62) 

However,  from  (54)  we  have  for  small  g  (i.e.,  g  <  g^) 

dr,  B3  -  £2/3  3^  B3  -  7r[Gi(-r;)  +  iAi(-r?)]  ,  (63) 

from  which  we  conclude  that  the  proper  choice  of  the  constant  C  is 
-Eq  and 

^(1)  -  -Eo  dg.  B3  .  (64) 

Having  used  the  behavior  of  B3  near  the  origin  to  find  the  constant  of 
proportionality  between  B3'  and  we  can  now  relate  the  field  strength 

Eq  to  the  amplitude  of  the  incident  whistler  wave.  Since  the  electric 
potential  of  the  incident  whistler  wave  is  represented,  in  the  mode 
conversion  solution,  by  B3  for  \g\-*  arg(?)-7r,  we  use  the  asymptotic  form 


of  and  (64)  to  make  this  connection.  Asymptotically  3^  -  dg.  B^CO) 

and  from  (40) 

dg  B^(0)  ~  i^2  r(l-a)  ei’fa  „-(l-a)ew/2  ^  (65) 

so  that 

-Eo  dg  B3  -  r(l/2+i^o)  Eo  w-(l--a)ew/2  .  (66) 

Since,  as  discussed  in  Sec.  IV,  W (^■^)exp(w/2)  represents  the  functional 
form  of  the  asymptotic  whistler,  the  coefficient  of  this  term  on  the  right 
hand  side  of  (66)  can  be  equated  to  the  electric  field  amplitude  of  the 
whistler  wave,  which  we  denote  as  ^2^'  Using  |  r(l/2+i^Q)|^  - 
Jn/cosh(n0Q) ,  we  then  obtain 

I  Eg  I  -  e-"'^°  [cosh(ff^o)]^^^  E^  _  (67) 

Jit 


Eq.  (67)  relates  the  asymptotic  value  of  the  whistler  electric  field  to 
the  value  of  the  driver  electric  field  of  Eq.  (60).  In  terms  of  physical 
quantities  the  parameter  has  the  form 


Po  - 


02 


2«(n2-w2)l/2 


k|L  . 


(68) 


VII.  Cone lus ions 


This  analytic  study  has  identified  a  new  mode  conversion  process  that 
can  occur  in  magnetized  plasmas  having  a  longitudinal  density  gradient 
(i.e.,  VnQ  X  B  -  0) ,  as  may  apply  in  the  auroral  ionosphere  and  long 
mirror  machines.  It  is  demonstrated  that  an  electrostatic  whistler  mode 
with  frequency  w  less  than  the  electron  gyrofrequency  (i.e.,  w  <  0) 
transfers  all  of  its  energy  to  a  hot  Bohm-Gross  mode  at  that  point  in  the 
density  profile  where  co  -  Wpg.  The  excited  hot  plasma  mode  propagates  away 
from  plasma  resonance  in  the  direction  of  decreasing  density  and 
eventually  transfers  its  energy  to  background  electrons.  Thus  the 
mechanism  described  here  provides  a  conceptual  link  between  the  presence 
of  large  amplitude  whistler  wave  activity  and  the  generation  of  magnetic 
field  aligned  fast  electrons.  Although  no  effort  is  devoted  in  the 
present  study  to  an  evaluation  of  the  quantitative  consequences  expected 
for  specific  physical  situations,  this  is  a  worthwhile  project  that  requires 
a  separate  survey. 

Our  major  effort  has  aimed  at  providing  an  analytical  description 
of  the  mode  conversion  process  and  an  evaluation  of  the  structure  of 
the  electric  field  excited  at  plasma  resonance.  The  solutions  of  the 
differential  equation  describing  the  mode  conversion  process  are 
divided  into  two  categories.  One  category  of  solutions  has  short 
wavelengths  and  is  related  to  Bohm-Gross  waves.  These  solutions  do  not 
exist  in  the  cold  plasma  limit.  Expressions  useful  for  evaluating  the 
asymptotic  properties  of  these  solutions  are  given  in  (20. a, b)  while 
the  series  expression  given  in  (21)  is  useful  for  evaluating  the 
solutions  near  the  origin.  The  second  category  of  solutions  is 


asymptotically  represent  propagating  electrostatic  whistler  waves.  The 
analytic  expressions  given  in  (29)  and  (36)  in  terms  of  Kummer 
functions  are  useful  for  evaluating  these  solutions  away  from  the 
origin.  The  long  wavelength  solutions  are  evaluated  near  the  origin  by 
relating  them  to  the  series  expressions  of  the  short  wavelength 
solutions.  The  physical  process  of  mode  conversion  is  shown  to 
correspond  to  the  mathematical  procedure  of  analytic  continuation  which 
relates  a  long  wavelength  solution  defined  over  a  given  argument  range 
to  a  combination  of  short  and  long  wavelength  solutions  defined  over  a 
separate  argument  range . 

Although  the  details  of  the  analysis  presented  are  somewhat  involved, 
they  point  to  a  powerful  simplification  in  describing  the  mode  conversion 
process.  It  is  shown  rigorously  that  at  plasma  resonance  the  excitation 
of  the  Bohm-Gross  mode  results  as  if  it  were  driven  by  a  uniform  capacitor 
plate  electric  field,  Eq ,  as  has  been  invoked  previously  in  other  resonant 
excitation  problems  (Morales  and  Lee,  1974  ;  Shoucri  and  Kuehl,  1980). 

The  analytic  expression  given  in  Eq.  (67)  relates  the  effective  pump  field 
Eq  to  the  asymptotic  whistler  amplitude  E^,  providing  the  connection 
needed  to  apply  the  capacitor  plate  model  to  an  incident  wave  problem. 

The  plasma  properties  appear  through  the  single  parameter  “ 
k^Lfi2/[2w(n2-w2)i/2] . 

Finally,  we  note  that  the  analysis  presented  here  complements  and 
extends  a  previous  study  (Maggs,  et.al.,  1984)  devoted  to  the  properties 
of  resonant  excitation  of  a  plasma  driven  externally  at  frequencies  larger 
than  the  electron  gyrofrequency  (i.e.,  w  >  fJ)  . 
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Appendix  I : 


In  this  appendix  a  series  representation  of  the  solutions  Aj  is 


derived  by  expanding  the  integrand  of  (17)  in  a  series  in  o^/a  and 


integrating  term  by  term.  Writing 


(I.l) 


where  the  coefficients  c  are 

m 

00 


c  -  >  (-)"  f  f  “+  ]  . 

n-0 


(1.2) 


and  f  “  1  “  — 2l! -  is  the  familiar  binomial  coefficient,  gives 

^  n! (a-n) ! 


A.  -  p .  £ 

J  J 


-2(o-l)/3 


^  c^  Oq™  g(f7,Q-m)  ,  (1.3) 


where 


g(r?,Q-m) 


a-m-2 


exp[  -(c7®/3  +  or?)]  .  (1.4) 


C(Aj) 


The  factors  p^  in  (1.3)  are  the  same  as  defined  for  the  as3nnptotic 


representation  of  the  functions  A^  given  in  (20. a, b).  Expanding 


exp  (-at?)  in  a  series  then  gives 


gin, a)  -  )  (-)^  I(Q!,k)  , 

k  k! 


(1.5) 


where  the  functions  g(r7,a)  are  the  same  as  used  by  Rabenstein  (1958),  and 


I(a,k)  “  dcr  o'*  ^  ^  exp  (-cr®/3) 


(1.6) 


C(Aj) 


The  integral  in  (1.6)  can  be  evaluated  by  making  the  change  of  variable, 


y  -=  0^/3.  The  contour  corresponding  to  then  transforms  into  one  which 
starts  at  infinity  with  argument  2ir  +  27r(2-j),  encircles  the  origin  in  the 
counter-clockwise  direction,  and  ends  at  infinity  with  argument  4?r  + 
27r(2-j).  Such  contours  correspond  to  Hankel's  contour  (Erdelyi,  et.al., 
1953)  so  that  (1.6)  can  be  written  in  the  form, 

I^(a,k)  -  27ri  exp(,ri(l-p)(7-2j))  /  (3^  r(p)]  ,  (1.7) 

where  p  -  (4-a-k)/3.  Using  (1.7),  (1.5)  and  (1.3)  then  gives  a  series 
expression  for  valid  in  the  finite  ?--plane  for  all  values  of  arg(?-) , 


A.  -  £ 
J 


-2(a-l)/3 


Pj  )l)l 


-)^  I,(tt-m,k)  /.  (1.8) 
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Appendix  II 


In  this  appendix  the  functional  form  of  the  electric  field  of  the  mode 


conversion  solution  is  obtained  in  terms  of  the  Airy  function,  Ai ,  and 


the  driven-Airy  function,  Gi.  Starting  with  equation  (53)  which  gives 


in  terms  of  the  derivatives  of  the  functions  A ^ '  with  respect  to  a, 


namely. 


B3'  -  +  aa(e-i2»ra.A^. )  +  _ 


(II. 1) 


and  using  the  m-0  term  of  the  series  representation  of  the  functions  A^ 


given  by  Eq.  (21)  or  (1.8),  we  obtain 


•  -  [a^ln  p  -  2  In  £]A  '  +  p.  V  8^1  ryk-l  , 

^  -^3  ^  ^  k^l  (k-D!  ^ 


(II. 2) 


where 


7ri(7-2j)  +  In  3  +  1^  , 


(II. 3) 


with  V’(p)  “  d[ln  r(p)]/dp  and  where  1^  is  given  by  (22)  or  (1.7). 


Defining  the  coefficients 


n  00 

I.  -  _ 1)  (-)^  V’(p)i.(p) 

J  3  k-1  (k-D!  J 


(II. 4) 


we  can  then  obtain  from  (II. 1)  using  (51)  and  (II. 2) 


B^'  -  [2e-i2»ra.A^.  4-A2')/3  +  [C^  +  +  C^]  .  (II. 5) 


Since  Ij  is  inversely  proportional  to  r(^)  ,  where  p  -  (3-k)/3  when  q-1 , 


the  factors  Cj  contain  the  term 


T  -  v>(p)/r(p)  -  r'ip)/rHp)  ■ 


(II. 6) 


To  evaluate  this  term  it  is  convenient  to  split  the  sum  over  k  into  three 


separate  sums  where 


in  case  1 : 


k  -  3£  ;  fc  -  0,1,2, . . . 


(II. 7. a) 


in  case  2; 


k-3n+l  ;  n-1,2,3,...  , 


(II. 7. b) 
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can  be  written  from  (II.  5)  as 

“  f(-r?)[  2nx  (2e'’''^/^  +  +  00^^(273))  +  cc^fi(-r7) 


■§(■'?)[  2?ri  c„  (2e’'’^/^  +  e"’^^/^)  +  CC^  \i(l/3)]  +  cc  fi(-fj) 


-  b(-fj)  , 


(11.18) 


where  c  -  1  +  .  Notii.g  that  c  -  0,  then  gives 


B3'  “  27ri  (l+e'^’f/3)c^f  (-r?)  -  {\+e^'^/^)c^g{-r})  ^  -  h(-r7)  ,  (11.19) 


From  the  definition  of  the  Airy  functions  (Antosiewicz ,  1970) 


Ai(x)  -  C3^f(x)  -  C2g(x)  , 


(II. 20. a) 


Gi(x)  -  3'i''2  +  C2g(x)]  +  x’ih(x)  . 


(II.20.b) 


B3'  can  be  simply  written  using  (11.19)  and  (II. 20. a, b)  as 


B3'  -  7r[  Gi(-fj)  +  iAi(-r?)]  . 


(11.21) 


FIGURES 


Figure  1.  An  electrostatic  whistler  with  frequency  w,  less  than  the 
electron  gyrofrequency  n,  propagates  to  plasma  resonance  (i.e.,  the  point 
in  the  plasma  where  tj-ajpg)  located  at  ^=0 .  The  gradient  in  the  plasma 
density  is  co-linear  with  the  magnetic  field. 


Figure  2.  The  integration  plane  is  cut  by  branch  lines  extending  from 
the  branch  points  s-=±i0^  to  infinity  parallel  to  the  real  axis.  The 
contours  defining  the  long  wavelength  solutions  to  (6)  all  start  at  the 
point  s  -  -^^2-  Examples  of  the  contours  for  this  class  of  solutions 
which  are  related  to  the  solutions  of  the  reduced  equation  (8)  are 
,  and  contours  are  shown  as  solid  lines  when  they  lie  in  the 

principal  Riemann  sheet  and  dashed  lines  when  they  lie  in  an  adjacent 
sheet.  The  contour  labeled  A^,  which  extends  to  infinity  in  sectors  2 
and  3,  is  an  example  of  a  short  wavelength  solution. 

Figure  3.  The  argument  range  of  f-  over  which  the  asymptotic 
expressions  for  Aj  given  in  (20. a)  and  (20. b)  are  valid  are  shown  for 
each  saddle  point  a  -  ±  The  function  behavior  at  values  of 

arg(?)  integer  multiples  of  ?r/3  is  also  indicated.  Specifically,  A^ 
behaves  isymptotically  as  a  left  to  right  propagating  wave  at 
arg(?-)  -  0. 


Figure  4.  The  principal  argument  ranges  for  the  functions  B^  (j  -  1,2,3  ) 
are  shown  in  relation  to  those  for  A^  .  Asymptotically,  the  functions 
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and  represent  left  to  right  propagating,  long  wavelength  waves  at 
arg(?')  -  0  and  tt  ,  respectively.  The  function  asymptotically  contains 
waves  propagating  in  both  directions  for  arg(^)  -  -tt.  The  argument  range 
of  the  function  B^  is  not  restricted.  At  arg(^)  -  0 ,  B^  contains  both 
left  and  right  propagating  waves. 

Figure  5.  a).  The  defining  contour  for  the  function  B^  can  be  distorted 

to  coincide  with  the  contours  defining  B^,  and  A^ .  Cauchy's  Theorem 

can  then  be  used  to  produce  a  relation  among  these  solutions  useful 
for  analytically  continuing  . 

b) .  A^  can  be  related  to  and  a  combination  of  B^  and  B^  by 
distorting  the  defining  contour  and  applying  Cauchy's  Theorem.  In 
addition  to  pulling  the  contour  in  towards  the  point  s  -  '^^2  ’  lower 

half  of  the  contour  for  is  pulled  around  the  point  s  -  in  the 

counter-clockwise  direction  and  back  towards  s  -  ‘^^2  passing  onto 

the  adjacent  Riemann  sheet. 
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